************************************************************************
************************************************************************
****             German Colonial State-Building
************************************************************************
************************************************************************

clear
cd "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/"

set more off



******************************
*** 1. Load Data
******************************

use "grid_panel_replication.dta"



******************************
*** 2. Main Analysis
******************************

* global controls
global controls "grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit##c.price_rubber_l i.year"

* Table 1
fvset base 1900 year 
* 5.1.1 DV: Stations, Station Type 
eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/main_stations.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear

* Output for Figure 2, need to run R code
logit station_count_nomiss $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
margins, at(battle_index_l=(-2.3 (0.1) 2.2)) post
parmest,saving(effect_battle)

logit station_count_nomiss $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
margins, at(stratval_n_correct_l=(-0.655 (0.1) 5.3419)) post
parmest,saving(effect_strat)


** Table 2
margins, at(battle_index_l=(.1145866 (.5232492) 2.2))
margins, at(stratval_n_correct_l=(-.2766822 (1.33) 5.3419))
margins, at(road_length=(0 (.2892523) 1.159514))
margins, at(ethnic_askari_number=(4 (92.99544) 315))
margins, at(dist_border=(91.78415 (100.7937) 394.9866))
margins, at(pop_density=(1.5 (.7561302) 4))







******************************
*** 3. Appendix
******************************

* Table 1 Summary Stats
estpost summarize grid_area cara_leng road_length centr_tribe1 pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast station_count_nomiss  near_station_dist battle_count_nomiss near_battle_dist battle_index soil_quality price_cotton rubber_suit price_rubber stratval_n_correct
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/summary_colonial.tex", cells("count(fmt(2)) mean(fmt(2)) sd(fmt(2)) min(fmt(2)) max(fmt(2))") noobs nonumber label title("Summary statistics") replace
eststo clear 


* Table 2: station distance with spatial lag
fvset base 1900 year
eststo:reg near_station_dist near_station_dist_splag_l $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:reg near_station_dist near_station_dist_splag_l $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:reg near_station_dist near_station_dist_splag_l $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_distance.tex",se scalars(r2 ll aic) label title(Distance to Nearest Station, OLS) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 3: Alternative Station Type
fvset base 1900 year
eststo:oprobit stat_type_code_nomiss_4 $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_4 $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_4 $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stationtype_other.tex",se scalars(ll aic) label title(Station Type, Ordered Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear



* Table 8: main models, ethnic centralization
eststo:logit station_count_nomiss $controls centr_tribe1 battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls centr_tribe1 battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls centr_tribe1 battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls centr_tribe1 battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls centr_tribe1 battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls centr_tribe1 battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_ethnic_centralization.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 9: main models, ethnic centralization +interaction 1
eststo:logit station_count_nomiss $controls c.centr_tribe1##c.battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls c.centr_tribe1##c.battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls c.centr_tribe1##c.battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls c.centr_tribe1##c.battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls c.centr_tribe1##c.battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls c.centr_tribe1##c.battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_ethnic_centralization_int1.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 10: main models, ethnic centralization + interaction 2
eststo:logit station_count_nomiss $controls c.centr_tribe1##c.stratval_n_correct_l battle_index_l ,r cluster(pu_id)
eststo:logit station_count_nomiss $controls c.centr_tribe1##c.stratval_n_correct_l2 battle_index_l2 ,r cluster(pu_id)
eststo:logit station_count_nomiss $controls c.centr_tribe1##c.stratval_n_correct_l5 battle_index_l5 ,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls c.centr_tribe1##c.stratval_n_correct_l battle_index_l ,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls c.centr_tribe1##c.stratval_n_correct_l2 battle_index_l2 ,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls c.centr_tribe1##c.stratval_n_correct_l5 battle_index_l5 ,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_ethnic_centralization_int2.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 11: Controlling for missions 
fvset base 1900 year
eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l mission_count,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2 mission_count,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5 mission_count,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/main_stations_missions.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear



* Table 12: non-endogenous variables
fvset base 1900 year 
* 5.1.1 DV: Stations, Station Type 
eststo:logit station_count_nomiss grid_area cara_len disease_count  elevmean annual_precipitation annual_mean_temp dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit##c.price_rubber_l i.year battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss grid_area cara_leng disease_count elevmean annual_precipitation annual_mean_temp dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit##c.price_rubber_l i.year battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss grid_area cara_leng disease_count elevmean annual_precipitation annual_mean_temp dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit##c.price_rubber_l i.year battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/main_stations_nonendo.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 13: battle index sub-components

fvset base 1900 year 
eststo:logit station_count_nomiss $controls battle_count_nomiss_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_count_nomiss_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_count_nomiss_l5 stratval_n_correct_l5,r cluster(pu_id)
eststo:logit station_count_nomiss $controls near_battle_dist_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls near_battle_dist_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls near_battle_dist_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_battle_count.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear



* Table 14: main models, ethnic region FE

set matsize 1000

eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l i.region,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2 i.region,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5 i.region,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l stratval_n_correct_l i.region,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l2 stratval_n_correct_l2 i.region,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l5 stratval_n_correct_l5 i.region,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_region_fe.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear

* Table 15: main models, ethnic group FE
eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l i.ethnicfe,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2 i.ethnicfe,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5 i.ethnicfe,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l stratval_n_correct_l i.ethnicfe,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l2 stratval_n_correct_l2 i.ethnicfe,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 $controls battle_index_l5 stratval_n_correct_l5 i.ethnicfe,r cluster(pu_id)

set matsize 800

esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_ethnic_fe.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 16: Grid Fixed Effects
xtset pu_id
eststo:xtreg station_count_nomiss $controls battle_index_l stratval_n_correct_l i.year,fe r cluster(pu_id)
eststo:xtreg station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2 i.year,fe r cluster(pu_id)
eststo:xtreg station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5 i.year ,fe r cluster(pu_id)
eststo:xtreg stat_type_code_nomiss_1 $controls battle_index_l stratval_n_correct_l i.year ,fe r cluster(pu_id)
eststo:xtreg stat_type_code_nomiss_1 $controls battle_index_l2 stratval_n_correct_l2 i.year,fe r cluster(pu_id)
eststo:xtreg stat_type_code_nomiss_1 $controls battle_index_l5 stratval_n_correct_l5 i.year,fe r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_full_fe_year.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 17: Correlated Random Effects Models
xtset pu_id
eststo:xtreg station_count_nomiss $controls battle_index_l_bw battle_index_l_wi stratval_n_correct_l_bw stratval_n_correct_l_wi,re r cluster(pu_id)
eststo:xtreg station_count_nomiss $controls battle_index_l2_bw battle_index_l2_wi stratval_n_correct_l2_bw stratval_n_correct_l2_wi,re r cluster(pu_id)
eststo:xtreg station_count_nomiss $controls battle_index_l5_bw battle_index_l5_wi stratval_n_correct_l5_bw stratval_n_correct_l5_wi,re r cluster(pu_id)
eststo:xtreg stat_type_code_nomiss_1 $controls battle_index_l_bw battle_index_l_wi stratval_n_correct_l_bw stratval_n_correct_l_wi,re r cluster(pu_id)
eststo:xtreg stat_type_code_nomiss_1 $controls battle_index_l2_bw battle_index_l2_wi stratval_n_correct_l2_bw stratval_n_correct_l2_wi,re r cluster(pu_id)
eststo:xtreg stat_type_code_nomiss_1 $controls battle_index_l5_bw battle_index_l5_wi stratval_n_correct_l5_bw stratval_n_correct_l5_wi,re r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_main_re.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 18: Lagged DV

set matsize 1000
fvset base 1900 year
eststo:oprobit stat_type_code_nomiss_1 stat_type_code_nomiss_l $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 stat_type_code_nomiss_l $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:oprobit stat_type_code_nomiss_1 stat_type_code_nomiss_l $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/stations_ldv_fe.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear

* Table 19: Instrumental Variables
tab year, gen(year_dum)
eststo:ivreg2 station_count_nomiss year_dum* grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast soil_quality price_cotton_l soil_price_int rubber_suit price_rubber_l rubber_price_int stratval_n_correct_l (battle_index_l = group_w_no ethn_viol ethn_def),r cluster(pu_id) 
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/iv_models.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear




* Table 20: Plantations
* collapse data to create cross-section
collapse rubber_count cotton_count station_count_nomiss grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast soil_quality price_cotton rubber_suit price_rubber battle_index stratval_n_correct loglight10, by(pu_id)

eststo:nbreg rubber_count grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast rubber_suit station_count_nomiss,r cluster(pu_id)
eststo:nbreg cotton_count grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast soil_quality station_count_nomiss,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/plantations.tex",se scalars(ll aic) label title(Number of Plantations, Logit) mtitles("Rubber" "Cotton") star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 21: Nightlights
eststo:reg loglight10 grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast soil_quality rubber_suit,r cluster(pu_id)
eststo:reg loglight10 grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean annual_precipitation annual_mean_temp dist_border dist_coast soil_quality rubber_suit station_count_nomiss,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/longterm_lights.tex",se scalars(ll aic) label title(Number of Plantations, Logit) mtitles("Rubber" "Cotton") star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear






* Table 4: Alternative Grid Cell Models
use "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/grid_panel_50_1.dta"

* global controls
global controls "grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean precip_mean temp_mean dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit_y##c.price_rubber_l i.year"

* Main Results
fvset base 1900 year 
* 5.1.1 DV: Stations, Station Type 
eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/main_stations_alt1.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear


* Table 5: Alternative Grid Cell Models 2
use "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/grid_panel_50_2.dta"

* global controls
global controls "grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean precip_mean temp_mean dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit_y##c.price_rubber_l i.year"

* Main Results
fvset base 1900 year 
* 5.1.1 DV: Stations, Station Type 
eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/main_stations_alt2.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear



* Table 6: Alternative Grid Cell Models 40km
use "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/grid_panel_40.dta"

* global controls
global controls "grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean precip_mean temp_mean dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit_y##c.price_rubber_l i.year"

* Main Results
fvset base 1900 year 
* 5.1.1 DV: Stations, Station Type 
eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/main_stations_alt40.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear





* Table 7: Alternative Grid Cell 60km
use "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/grid_panel_60.dta"

* global controls
global controls "grid_area cara_leng road_length disease_count ethnic_askari_number ethnic_number pop_density elevmean precip_mean temp_mean dist_border dist_coast c.soil_quality##c.price_cotton_l c.rubber_suit_y##c.price_rubber_l i.year"

* Main Results
fvset base 1900 year 
* 5.1.1 DV: Stations, Station Type 
eststo:logit station_count_nomiss $controls battle_index_l stratval_n_correct_l,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l2 stratval_n_correct_l2,r cluster(pu_id)
eststo:logit station_count_nomiss $controls battle_index_l5 stratval_n_correct_l5,r cluster(pu_id)
esttab using "~/Dropbox/German_Colonialism/Statebuilding_Paper/Replication_Archive/Tables/main_stations_alt60.tex",se scalars(ll aic) label title(Station Presence, Logit) mtitles() star(+ 0.10 * 0.05 ** 0.01 *** 0.001) replace
eststo clear






















